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(-H \ Using both cosmic microwave background and large scale structure data, we put stringent 

bounds on the possible cold dark matter, baryon and neutrino isocurvature contributions to 
primordial fluctations in the Universe. Neglecting the possible effects of spatial curvature, 
tensor perturbations and reionization, we perform a Bayesian likelihood analysis with nine 
free parameters, and find that the relative amplitude of the isocurvature component cannot 
be larger than about a few % for uncorrelated models. On the other hand, for correlated 
adiabatic and isocurvature components, the ratio could be slightly larger. However, the cross- 
correlation coefficient is strongly constrained, and maximally correlated/anticorrelated models 
are disfavored. 
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Thanks to the tremendous developments in observational cosmology during the last few 
years, it is possible to speak today of a Standard Model of Cosmology, whose parameters are 
known within systematic errors of just a few percent. Moreover, the recent measurements of 
both temperature and polarization anisotropies in the cosmic microwave background (CMB) has 
opened the possibility to test not only the basic paradigm for the origin of structure, namely 
inflation, but also the precise nature of the primordial fluctuations that gave rise to the CMB 
anisotropies and the density perturbations responsible for the large scale structure (LSS) of the 
Universe. 

The simplest realizations of the inflationary paradigm predict an approximately scale in- 
variant spectrum of adiabatic and Gaussian curvature fluctuations, whose amplitude remains 
constant outside the horizon, and therefore allows cosmologists to probe the physics of infla- 
tion through observations of the CMB anisotropies and the LSS matter distribution. However, 
this is certainly not the only possibility. Multiple-field inflation predicts that, tog ether with 
the adiabatic component, there should also be entropy or isocurvature perturbations* * *****, 
associated with fluctuations in number density between different components of the plasma be- 
fore decoupling, with a possible statistical correlation between the adiabatic and isocurvature 
modes ' . Baryon and cold dark matter (CDM) isocurvature perturbations were proposed long 



agolHlas an alternative to adiabatic perturbations. Recently, two other modes, neutrino isocur- 
vature density and velocity perturbations, have been added to the listr*! Topological defects also 
predict a significant isocurvature component. Moreover, it i s well known that entropy pertur- 
bations seed curvature perturbations outside the horizonl^^, so it is possible that a significant 
component of the observed adiabatic mode could be stron gly correlated with an isocurvature 
mode. Such models are generically called curvaton models'^, and are now widely studied as 
an alternative to the standard paradigm. Furthermore, isocurvature modes typically induce 
non-Gaussian signatures in the spectrum of primordial perturbations. 

I will review here our recent results'*^ constraining the various isocurvature components, 
using data from the temperature power spect rum and temperature-polarization cross-correlation 
recently measured bv the WMAP satellite ^^ from the small-scale temperature anisotropy 
probed by ACBAR^"", and from t he m atter power spectrum measured by the 2-degree-Field 
Galaxy Redshift Survey (2dFGRS)^l We do not use the data from Lyman-a forests, since 
they are based on non-linear simulations carried under the assumption of adiabaticity. We will 
not assume any specific model of inflation, or any particular mechanism to generate the pertur- 
bations (late decays, phase transitions, cosmic defects, etc.), and thus will allow all five modes 
- adiabatic (AD), baryon isocurvature (BI), CDM isocurvature (CDI), neutrino isocurvature 
density (NID) and neutrino isocurvature velocity (NIV) - to be correlated (or not) among ea- 
chother, and to have arbitrary tilts. However, we will only consider the mixing of the adiabatic 
mode and one of the isocurvature modes at a time. This choice has the advantage of restricting 
considerably the parameter space, and reflects the fact that most of the proposed mechanisms 
for the generation of isocurvature perturbations lead to onl y one mode. The first bounds on 
isocurvature perturbations assumed uncorrelated modes'^, but recently also correlated ones 
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were considered in Refs. 1 ' ' ' '. 

For simplicity, we also neglect the possible effects of spatial curvature, tensor perturbations 
and reionization. Therefore, each model is described by nine parameters: the cosmological 
constant Q\, the baryon density ujb = Q~Bh 2 , the cold dark matter density w C( j m = f2 c d m /i 2 , the 
overall normalisation A, the relative amplitude a of the isocurvature to adiabatic modes, and 
their correlation /3, the adiabatic and isocurvature tilts, n at j and n; so , and finally a matter bias b 
associated to the 2dF power spectrum. We generate a 5-dimensionnal grid of models (^4, a, /3, 
and b are not discretized) and perform of Bayesian analysis in the full 9-dimensional parameter 
space. At each grid point, we store some C\ values in the range < I < 1800 and some P(k) 
values in the range probed by the 2dF data. The likelihood of each model is then computed 
using the software or the detailed information provided on the experimental websites, using 1398 
points from WMAP, 11 points from ACBAR and 32 points from the 2dFGRS. For parameter 
values which do not coincide with grid points, our code first performs a cubic interpolation of 
each power spectrum, accurate to better than one percent, and then computes the likelihood 
of the corresponding model. From the limitation of our grid, we impose a flat prior on the 
isocurvature tilt: rii so > 0.6 - inflation predicts the tilts of both adiabatic and isocurvature 
modes to be close to one. The other grid ranges are wide enough in order not to affect our 
results. 

For t he theoretical analysis, we will use the notation and some of the approximations of 
Refs.™4 During inflation more than one scalar field could evolve sufficiently slowly that their 
quantum fluctuations perturb the metric on scales larger than the Hubble scale during inflation. 
These perturbations will later give rise to one adiabatic mode and several isocurvature modes. 
We will restrict ourselves here to the situation where there are only two fields, <fii and </>2, 
and thus only one isocurvature and one adiabatic mode. The evolution during inflation will 
draw a trajectory in field space. Perturbations along the trajectory (i.e. in the number of e- 
folds N) will give rise to curvature perturbations 1Zk = 5Nk = H 5k = H5pk/p on comoving 
hypersurfaces, while perturbations orthogonal to the trajectory will give rise to gauge invariant 



entropy (isocurvature) perturbations, S = 5\n{rii/rij) = 5pi/{pi + pi) — 5pj/(pj +Pj)? 

In order to relate these perturbations during inflation with those produced during the ra- 
diation era, one has to follow the evolution across reheating. During inflation we can always 
perform an instantaneous rotation along the field trajectory and relate perturbations in the 
fields, 5<pi and 5(p2, with perturbations along and orthogonal to the trajectory, 5a and 5s, 
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In this cas e, t he curvature and entropy perturbations can be written, in the slow-roll approxi- 





mation, as 1 
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The problem however is that, contrary to the case of purely adiabatic perturbations, the 
amplitude of the curvature and entropy perturbations does not remain contant on superhorizon 
scales, but evolve with time. In particular, due to the conservation on the energy momentum 
tensor, the entropy perturbation seeds the curvature perturbation, and thus their amplitude 
during the radiation dominated era evolves according to 



K k = a{t)HS k 
Sk = f3(t)HS k , 



(4) 
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where a and (3 are time-dependent functions characterizing the evolution during inflation and 
radiation eras. A formal solution can be found in terms of a transfer matrix, relating the 
amplitude at horizon crossing during inflation (*) with that at a later time during radiation, 
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where the transfer functions are given by 



l3(t')H(t')dt' 



Tss(t,U) = exp 

Tjis(t,U)= f a(t')H(t')T ss (t',U)dt'. 
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Note that in the absence of primordial isocurvature perturbation, 5* = 0, the curvature per- 
turbation remains constant and no isocurvature perturbation is generated during the evolution. 
This is the reason for the entries T-pm = 1 and Tg-ji = 0, respectively, in the transfer matrix. 

Since <pi and fa are essentially massless during inflation, we can treat them as free fields 
whose fluctuations at horizon crossing have an amplitude Sfa = {H k /yj2k^) ej(k), where H k is 
the rate of expansion at the time the perturbation crossed the horizon (k* = aH), and ej(k) are 
gaussian random fields with zero mean, (ej(k)) = and (ej(k)e^(k')) = 5ij <5(k — k'). Now, since 
5a k and 5s k and just rotations of the field fluctuations, they are also gaussian random fields of 

a For instance, entropy perturbations in cold dark matter during the radiation era can be computed as <Scdi = 
5 m — 3<5 r /4. 



amplitude Hk/V2k' i . However, the time evolution Q will mix those modes and will generically 
induce correlations and non-gaussianities. 

Therefore, the two-point correlation function or power spectra of both adiabatic and isocur- 
vature perturbations, as well as their cross-correlation, can be parametrized with three power 
laws, i.e. three amplitudes and two spectral indices, 

t.3 / h \ "ad- 1 

*)4( 5! »^ , S • (9) 

A^(fc) = ^(^5)=^i?cosA (^-J 

where fcg is some pivot scale, typically of the order of the present horizon. Note that we have 
assumed here that the correlation coefficient cos A is scale-independent. In general there will be 
a small scale dependence coming from that of t*(k) in the transfer functions l|7jl. We will ignore 
it for the moment. 6 

The CMB temperature fluctuations on large scales, the Sachs- Wolfe contribution, can be 
written as 

^ = -<!>£, + 2<Df a ° d = l -K r&A - ~«S rad , (10) 

in terms of the initial conditions during the radiation era, ^ad = (3/5)7£ ra d and 3>J.ad = 
— (l/5)<S ra d. However, in order to evaluate the temperature and polarization anisotropies seen in 
the CMB today, one has to calculate the radiation transfer functions for adiabatic and isocurva- 
ture perturbations, gf d (k), gf°(k) and gf°"(k). The angular power spectrum for each individual 
mode is therefore 
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and the total angular power spectrum is computed as^^ 

Q = Cf + B 2 Cf° + IB cos A Cf ori , (11) 

where A = 1 has been marginalized over, and B is now the entropy to curvature perturbation 
ratio during radiation era. , B = S raA /lZ vax ±. We will use here a slightly different notation, used 
before by other groups^ ' \ where A 2 = (1 — a) and B 2 = a, 



Q = (1 - a) Cf + a Cf ° + 2(3 yja(l - a) Q corr . (12) 

Here the parameter a runs from purely adiabatic (a = 0) to purely isocurvature (a = 1), 
while (3 defines the correlation coefficient, with (3 = +1/ — 1 corresponding to maximally corre- 
lated/anticorrelated modes. There is an obvious relation between both parametrizations: 
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B 2 /(l + B 2 ), /? = cosA. (13) 



In a recent paper*-^* this assumption was relaxed, but produced results that are not very sensitive to it. 
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Figure 1: a) the likelihood function of the isocurvature fraction a, for the three different types of uncorrected 
isocurvature modes (i.e., with the prior f3 = 0); b) the 1- and 2-ct contours of a and the cross-correlated mode 
coefficient 2/3-y/a(l — a), for the CDM isocurvature mode: the small (red) contours are based on all the data, 
with one flat prior n\ so > 0.6, while the large (green) ones show the situation before WMAP, with an additional 
prior lob < 0.037; c) same as b) for NID; d) same as b) for NIV. 



This notation has the advantage that the full parameter space of (a, 2[3y/a(l — a)) is contained 
within a circle of radius 1/2. The North and South rims correspond to fully correlated {(3 = +1) 
and fully anticorrelated (/3 = —1) perturbations, with the equator corresponding to uncorrelated 
perturbations ((3 = 0). The East and West correspond to purely isocurvature and purely adia- 
batic perturbations, respectively. Any other point within the circle is an arbitrary admixture of 
adiabatic and isocurvature modes. 

In order to compute the theoretical prediction for the C\ coefficients of the temperature and 
polarization power spectra, as well as the matter spectra P(k), for all four different components, 
we have used a CMB code developed by one of us (A.R.) which coincides, within 1% errors, 
with the values provided by CMBFAST for AD and CDI modes, and also includes the neutrino 
isocurvature modes as well as the cross-correlated power spectra. Note that the code defines the 
tilt of each power spectrum with respect to a pivot scale fco corresponding the present value of 
the Hubble radius, while CMBFAST uses k^ = 0.05 Mpc -1 . Thus, the comparison of our results 
with those of WMAP for the CDI mixed model is not straightforward. 
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Figure 2: Some best fit models compared to data. The left plots correspond to correlated CDM isocurvaturc 
modes, the right plots to neutrino density isocurvature modes. The parameters of these models can be seen in 
Table 1. From top to bottom, we show the C ; , C l and P (k) power spectra, as well as the contribution of 
each component: adiabatic (ad), isocurvature (iso), cross-correlated (corr). The last two components have been 
rescaled by a factor indicated in each figure. We also show the data points that we use throughout the analysis, 
from WMAP, ACBAR and 2dF. In the case of the matter power spectrum, we plot in addition the theoretical 
spectrum convolved with the experimental window function (conv), in order to allow for a visual comparison with 
the data (which is not possible directly with the unconvolved power spectrum, due to the shape of the window 
functions). Note that these two best-fit models are extremely close to the pure adiabatic one: the differences in 
the respective power spectra would be indistinguishable by eye. 



Table 1: The best fit values of cosmological parameters for both adiabatic, mixed-isocurvature and correlated 

models, with their y 2 P er degree of freedom. 



model 
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AD 


0.021 


0.12 


0.70 


0.95 


- 


- 


- 


1478.8/1435 


CDI 


0.023 


0.12 


0.73 


0.99 


1.02 


0.10 


- 


1478.3/1433 


BI 


0.023 


0.12 


0.73 


0.99 


1.02 


0.72 


- 


1478.3/1433 


NID 


0.023 


0.12 


0.73 


0.99 


0.95 


0.37 


- 


1478.2/1433 


NIV 


0.021 


0.12 


0.70 


0.95 


- 


0.0 


- 


1478.8/1433 


c-CDI 


0.022 


0.12 


0.71 


0.97 


1.23 


0.001 


1.0 


1478.0/1432 


c-BI 


0.022 


0.12 


0.71 


0.97 


1.23 


0.03 


1.0 


1478.0/1432 


c-NID 


0.022 


0.12 


0.73 


0.97 


1.04 


0.10 


0.26 


1477.7/1432 


c-NIV 


0.021 


0.12 


0.71 


0.95 


0.71 


0.03 


-1.0 


1477.5/1432 



Conclusions. Using the recent measurements of temperature and polarization anisotropies 
in the CMB by WMAP, as well as the matter power spectrum measured by 2dFGRS, one can 
obtain stringent bounds on the various possible isocurvature components in the primordial spec- 
trum of density and velocity fluctuations. We have considered both correlated and uncorrelated 
adiabatic and isocurvature modes, and find no significant improvement in the likelihood of a cos- 
mological model by the inclusion of an isocurvature component. For instance, a model allowed 
at the 2a level by the present observations (WMAP + ACBAR + 2dFGRS) with a correlated 
admixture of adiabatic and CDM isocurvature modes, with a = 0.46 and (3 = —0.3, and tilts 
n a d = 1.03 and n; so = 0.79, together with cosmological parameters $7a = 0.744, w c dm = 0.112 
and wb = 0.0246, has a x 2 per d.o.f. of 1.0349, while the corresponding c-CDI best fit model (see 
Table. 1) has x 2 l v = 1-0321, and the pure adiabatic best fit model has x 2 l v = 1-0305. There- 
fore, we conclude that for the moment there is no significant improvement in the likelihood of 
a model by the inclusion of a small admixture of isocurvature perturbations. In other words, 
the basic paradigm of single fiel d in flation passes for the moment all observational constraints 
with flying colors. It is expected'^ that in the near future, with better data, we will be able to 
constrain further a possible isocurvature fraction, or perhaps even discover it. 
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